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Abstract 

We check a recent proposal [H. Goto and K. Ichimura Phys. Rev. A 70, 012305 (2004)] for 
controlled phase gate through adiabatic passage under the influence of spontaneous emission and 
the cavity decay. We show a modification of above proposal could be used to generate the necessary 
conditional phase gates in the two-qubit Grover search. Conditioned on no photon leakage either 
from the atomic excited state or from the cavity mode during the gating period, we numerically 
analyze the success probability and the fidelity of the two-qubit conditional phase gate by adiabatic 
passage. The comparison made between our proposed gating scheme and a previous one shows 
that Goto and Ichimura's scheme is an alternative and feasible way in the optical cavity regime for 
two-qubit gates and could be generalised in principle to multi-qubit gates. 
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I. INTRODUCTION 



The stimulated Raman adiabatic passage (STIRAP) has been extensively studied in 
coherent population transfer By partially overlapping pulses in the counterintuitive 
sequence, we have the population efficiently transferred between two quantum states while 
almost not populating the intermediate level. Thus the effect of spontaneous emission from 
the intermediate level is negligible. Moreover, the STIRAP is independent of the pulse 
shape, which makes it robust against moderate fluctuations of experimental parameters Q • 
Because of the above mentioned merits, many schemes concerning quantum information 
processing (QIP) are based on the technique of STIRAP, such as single qubit rotation 
controlled NOT gate 4], controlled phase gate ^,3, arbitrary state controlled- unitary gate 
SWAP gate [8], generation of qubit entanglement [9] and so on. 

All the above mentioned QIP schemes are quite different from the dynamical ones flo| . 
which need a precise control of the Rabi frequency and the pulse duration. They are also 
quite different from the adiabatic geometric ones [HJ| , which depend on a controllable loop 
in the parameter space. In Ref. 5|], an attractive scheme to generate multi-qubit controlled 
phase gate by only three steps in an optical cavity has been presented. As in Refs. j^, Q, f|, 
it utilizes the cavity mode to couple different atomic qubits. Although the dark states used 
for adiabatic evolution include no components of atomic excited state, they usually contain 
one photon or multi-photon cavity state. That is to say, the dark states are decoherence-free 
states with respect to the spontaneous decay, but not immune to the cavity decay. However, 
if the laser intensity is much smaller than the atom-cavity coupling strength [7, Q], the 
population of the cavity state with non-zero photon can be very little. Nevertheless, with a 
weak laser intensity, the pulse duration time must be much prolonged in order to fulfill the 

n q 

adiabatic condition 12J , which would enhance the probability of dissipation. 

Therefore, it is natural for us to ask how well the quantum gating in Ref. 3 works 
under the influence of dissipation. In the present paper, we will show how to generate all 
the conditional phase gates in a two-qubit Grover search by the STIRAP technique, based 
on the main idea in Ref. To check the influence from spontaneous emission and the 
cavity decay, we will employ the quantum jump approach [12J, which uses wave function to 
describe the system subject to dissipation while master equation employs density matrix. 
Conditioned on no decay happening, the system's evolution is governed by a non-Hermitian 
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Hamiltonian. In some cases, quantum 
for systems subject to dissipation 



; ump approach could help us find analytical solutions 
131 ] . But in the present case, due to complexity, we 
will have to numerically simulate the gating process and investigate the dependence of the 
success probability and the gate fidelity on the spontaneous emission rate T and the cavity 
decay rate k. 

II. CONDITIONAL PHASE GATE WITHOUT DECAY 

The conditional phase gate is used for labeling target state in the Grover search algorithm 
Q, [isj], which means an addition of a it phase as a prefactor to the target state, but of 
nothing to other states. The target state labeling is a key step in Grover search. For items 
represented by the computational states \X) with X — 0, 1, N — 1, in a quantum register 
with n qubits, we have N = 2 n possible states. If the target state is |r), the conditional 
phase gate can be expressed as I T — I — 2|r)(r|, where / is the N x N identity matrix. By 
redefining the energy levels in Ref. |5|, we show that the second step in Ref. 5] to generate 
the two-qubit controlled phase gate is actually for a conditional phase gating to label the 
target state 1 0) 1 1 1 ) 2 - Atoms are fixed in the optical cavity, as shown in Fig. 1, for an atom j 
with five-level configuration, where levels \0)j and \a)j are coupled to the excited state \2)j 
by two lasers with Rabi frequencies and Q a j respectively, while the level |l)j is coupled 
to state \2)j by the cavity mode with the coupling constant gj. These three couplings are 
needed to construct a two-qubit or even a multi-qubit gate. The couplings of the levels \0)j, 
\a)j and |l)j to another excited state 1 3) j can be used to perform single qubit rotation as in 
Ref. 3. 

We encode the qubits in levels |0)j and |l)j, and for simplicity we focus our discussion on 
the case of two atoms, although our case is extendable to multi-atom cases. So our task is 
to accomplish I^^ju, i.e., adding a minus sign to the target states | « ) 1 1 j? ) 2 where i, j being 
0, 1, respectively and the subscripts are for different atoms. Consider the pulses regarding 
f2 0j i and fl a ^, the Hamiltonian is given by (assuming h = 1) 

H = f2 ,i|2)n(0| + ri (T) 2|2)22{c| + ^ ^ 9mO\ 2)mm ( 1 1 1 1 • (1) 

m=l,2 

where a is the annihilation operator for the cavity mode. As the cavity is initially 
in the vacuum state |0), |l)i|0) 2 |0) and |1) 1 |1) 2 |0) are unaffected by the Hamiltonian, 



while |0)i|0)2|0) and |0)i|l)2|0) are associated with the Hamiltonian's two dark states 
(i.e., eigenstates with zero eigenvalues) |A)o) oc 1 0) 1 1 0) 2 1 0) — Oo,i 1 1) 1 1 0) 2 1 1) » |A)i) oc 
S'i^ct^ |0> 1 1 1> 2 10) + 5'2^o,i|l)i| cr )2|0) - ^o,i^cr,2 1 1) 1 1 1) 2 1 1) respectively. The procedure con- 
sists of two STIRAP processes: (i) In the first STIRAP, fi CT)2 precedes O ,i; (h) The second 
STIRAP is a reverse process of the first one but with the phase regarding Q a ^ 2 added by n. 
n the adiabatic limit and on the condition that the Berry phase is equal to zero, we have 
5I (assuming f2 ,i, fi CTj2 having the same phase in (i) and g\ — # 2 — g) 

|0> 1 [0> 2 |0> [0> 1 |0> 2 |0> ^ |0> x |0> 2 [0> (2) 

|o>x|i>2|o> $ |i}ik) 2 |o) n — 10> 1 |1> 2 |0> 
|l)xi0} 2 |0) ^ I1> X |0> 2 |0> ^ |l}i|0} 2 |0) 
|l)i|l) 2 |0> ^ |l)i|l) 2 |0) ^ |l>x|l> 2 |0> 

Thus the conditional phase gate for labeling the target state |0)i|l) 2 is generated. All the 
other conditional phase gates can be generated by adding the NOT gate on both sides of 
the above two STIRAP processes. It is easy to see i|o)i|o) 2 = °'x,2^|o)i|i>2 <J x,2, ^|i)i|o) 2 = 
c r x,i cr x',2^|o>i|i)2 "x,2C r x,i ) J |i>i|i> 2 = o-x,i^|o)i|i)2 cr x,i, where cr Xti (i = 1,2) is the NOT gate acting 
on the ith atom, transforming states as |0)j — > — > |0)j. The NOT gate a x ^ can 
be obtained by coupling |0}j, \a)i and with the excited state |3)j, by three STIRAP 
processes: (1) Q aji precedes f^; (2) precedes Qo ti ; (3) fi ,j precedes Q^. The dark 
states associated with the above three steps are \D^) oc f2 CT) i|l)j— f2i,i|cr)j, \D^) oc ^i,i|0)j— 
Qoj\l)i, \D^) oc ilo^l^i — r2 CT .j|0)j respectively. In the adiabatic limit and with the condition 
of zero Berry phase, we obtain in the case of the two pulses in each STIRAP with tt phase 
difference, 



io), a | 0)i n id, a id,, 

ID, 3 Wh 3 !„>, n |o)„ 



(3) 



As a result, we realize the NOT gate. Alternatively, the NOT gate can also be reached by 
single qubit rotation as in Ref. [3J, which also needs 6 pulses. 
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III. CONDITIONAL PHASE GATE WITH DECAY 



We introduce the spontaneous emission and the cavity decay from now on and will nu- 
merically analyze their effects on the conditional phase gate by the quantum jump approach. 
As long as there is no photon leakage either from the atomic excited state or from the cavity 
mode during the gating period, the Hamiltonian in Eq. (1) becomes 



H cond = [fi ,i|2)ii(0|+a T , 2 |2) 22 H+ ga\2) mm {l\ + H.c.]-i-<Ja-i- ^ |2> mm (2|, (4) 

m=l,2 m=l,2 

where we have made g± — (? 2 — g. The above Hamiltonian can be written as a matrix in the 
subspace spanned by |0) 1 |1) 2 |0), |2) 1 |1) 2 |0), |l)i|l) 2 |l), |l)i|2) 2 |0), |1)i|<t) 2 |0>, |0)i|0) 2 |0), 
|2>i|0> 2 |0>, |1> 1 |0> 2 |1>, |l>i|0> 2 |0>, |l>i|l> 2 |0> 
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(5) 



Note that H cond is non-Hermitian, the norm of a state vector evolving under the corre- 
sponding Schrodinger equation decreases in general with time. For an arbitrary initial state 

Psuc{t) = mrWlnA^Ucond^t^U)) (6) 

defines the probability that no photon has been emitted at time t, where U conc i is the time 
evolution operator for H con d and the corresponding normalized state vector is 



\i/>(t)) = U cond (t,U)\ij(t t ))/y/P suc (t). 

The gate fidelity is 

5 



(7) 



F=|^(oo)|/| 0>1 | 1>2 |^))| 2 . (8) 

In our numerical simulation, we choose fi 01max = ff)2 , maj[ = 0.16g in order to reduce 
the population of single-photon cavity states during the gating process, i.e., |l)i|0) 2 |l) and 
1 1) i j 1)2! 1) - Similar to Refs. {3], 5], we suppose that the pulses have the Gaussian shape 
exp[-(t ± T/2 ± t ) 2 /2r 2 ] with T = 200/g, t = 30/p, and r = 40/g. In Fig. 2, we show 
the time evolution of the initial state §(|0)i|0) 2 + |0)i|l) 2 + |l>i|0) 2 + |l)i|l) 2 )|0). Ideally, 
the final state would be |(|0) 1 |0) 2 - |0)i|l) 2 + |l)i|0) 2 + |l)i|l) 2 )|0). By our numerical 
calculation, with k = O.lg and T = O.lg, the probability amplitudes for states |0)i|0) 2 |0), 
|0)i|l) 2 |0), |l)i|0) 2 |0), |l)i|l) 2 |0) are 0.4513, -0.4523, 0.5438, 0.5438 respectively, and the 
fidelity is 99.12% while the success probability is 84.55%. The reason for the amplitudes of 
the components |l)i|0) 2 |0) and |l)i|l) 2 |0) increased compared with those in the initial state 
is that they did not participate in the evolution, while the other two components dissipated 
in the evolution. So the relative weights of |l)i|0) 2 |0) and |l)i|l} 2 |0) are enlarged in the 
normalized final state. 

In order to see how the success probability P suc and the gate fidelity F depend on k and T, 
we plot P suc and F in each of Figs. 3, 4, 5, for the evolution from the initial states |0)i|l) 2 |0), 
|0)i|0) 2 |0), |(|0)i|0) 2 +|0) 1 |l) 2 +|l)i|0) 2 +|l) 1 |l) 2 )|0), respectively. From Fig. 3, we see that 
P suc is less affected by k than by T, while it is almost unaffected by V in Fig. 4. Since for 
the evolution of |0)i|l) 2 |0), the occupation probability of |l)i|l) 2 |l) is several times smaller 
than that of the excited atomic state, the spontaneous emission from the excited state |2) 
is more detrimental than cavity decay in Fig. 3. While for the evolution of |0)i|0) 2 |0) in 
Fig. 4, as the occupation in level |2)i is negligible, the cavity decay is the main detrimental 
effect in this case. For an arbitrary initial state, \i>{U)) = a|0)i|l) 2 |0) + /5|0)i|0) 2 |0) + 
7|l)i|0) 2 |0) +e|l)i|l) 2 |0), we have U cond (oo, UMiU)) = aV Psudl^oio) + /V-FUd^ooo) + 
7|l)i|0) 2 |0) + e|l)x|l) 2 |0), where l^oio) and \4>ooo) are, respectively, the final states after 
time evolution of the component states |0)i|l) 2 |0), |0)i|0) 2 |0), and P SU ci, P S uc2 are their 
corresponding success probabilities. Thus the success probability for any initial state is 
Psuc — W\ 2 P S uci + \/3\ 2 P S uc2 + |7| 2 + k| 2 - Usually both P suc i and P SUC 2 are smaller than 1. As 
a result, for the same T and k, the less the initial population in components |0)i|l) 2 |0) and 
|0)i|0) 2 |0), the higher the success probability for the gating process. Moreover, from Figs. 
3 and 4, we know that the fidelity is very high for the initial state |0)i|l) 2 |0) or |0)i|0) 2 |0), 
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which implies l^oio) — — 1 0) 1 1 1)2 1 0) , |^ooo) — 1 0) 1 1 0) 2 1 0) , while the relative weight ratio of 
the four component states |0)i|l) 2 |0), |0)i|0) 2 |0), |l)i|0) 2 |0), |l)i|l) 2 |0) is approximately 
{—o.y/P suc \, (3yJ Psuci-, 7, e) , deviated from the ideal case (—a, /3, 7, e). Obviously, the 
fidelity for the state \ip(ti)) would be very high in the following two cases: 1) 7 = e = 0, 
Psud — P S uc2] 2) 7, e ^ 0, P sucl ~ P suc2 ~ 1. Therefore, it is easy to understand in Fig. 5 
why P suc is increased in most regions while F is decreased by comparing with those in Figs. 
3 and 4 with the same T and k. 

We also did the simulation for the NOT gate (not shown), which shows that with intense 
short pulses and the vacuum cavity state, the fidelity and the success probability are both 
near unity even under large k and T. So the success probability and the fidelity for the 
other three conditional phase gates would have almost the same dependence on k and V as 
in above figures. 



IV. DISCUSSION AND CONCLUSION 



Combined with single qubit Hadamard gate, which is achievable by a series of single qubit 
rotations as in Ref. [3J, the performance of the two-qubit Grover search by the conditional 
phase gates is straightforward Jj], 15]. Our scheme has the following advantages compared 
with Ref. Atoms are fixed in the cavity, which is easier to control than the case with 
atoms going through the cavity with a certain velocity. Moreover, all the operations are 
based on STIRAP, so they are very robust to fluctuation of experimental parameters. 

We now briefly discuss the experimental feasibility. It has been experimentally reported 
that g 



34 x 2vr MHz, k = 4.1 x 2tt MHz and T = 2.6 x 2tt MHz [16j. If we adopt the 
above numbers, i.e., n/g ~ 0.12, Y/g ~ 0.077, the operation time is of the order of 10 -6 s, 
as can be seen from Fig. 2(a). In this case, for the initial states |0)i|l) 2 |0), |0)i|0) 2 |0) and 
|(|0)i|0) 2 + |0)i|l) 2 + |l) 1 |0) 2 + |l)i|l) 2 )|0), we may obtain that P suc ~ 73.1%, 63.8%, 84.2%, 
and F ~ 99.9%, 100%, 99.0%, respectively. The fidelities are high enough, while the success 
probabilities need to be improved by reducing kJq and Y/g in the future experiment. 



We have noticed a previous smart scheme 



171 ] for two-qubit phase gate, based on the 



STIRAP technique, by an almost same configuration as in the present paper. Assisted 



by Zeno effect, the state involved in adiabatic evolution in Ref. 



171 ] is immune to both 



the cavity decay and spontaneous emission. From their simulation for a STIRAP process 
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with fixed parameters k = O.lg and T = O.lg, we see the variation of the maximum Rabi 
frequency and the detuning of the corresponding transitions of the atomic levels can lead to 
the maximum success probability and the high fidelity for the population transfer jl8!]. In 
contrast, we did simulations in our work for successive two STIRAP processes without any 
change of the detuning and the pulse Rabi frequency ljj. Given all the facts above and 
by comparing the simulation results, we consider that the proposal presented in Ref. Qjjs 



not less effective than that in Ref. 



17] in the presence of dissipation. Moreover, Ref. 



is applicable to multi-qubit phase gates, which would be used in multi-qubit Grover search 
algorithm. But with more qubits involved, the dark states containing components with more 
photons, would be more sensitive to dissipation. However, with very high Q cavity, Ref. 
can be in principle feasible in the case of multi-qubit phase gates. 

In summary, we have shown a generation of four conditional phase gates required by a 
two-qubit Grover search algorithm by means of the STIRAP technique, in the presence of 
dissipation. To meet the requirement of adiabatic condition, the multiplication of the Rabi 
frequency and the implementation time should be much larger than unity. But the larger 
the Rabi frequency, the more the occupation probability in the cavity mode, yielding more 
probability of the cavity decay. On the other hand, the less Rabi frequency would require a 
longer implementational time, which would also enhance the probability of dissipation. In 
this sense, it is of importance to find suitable numbers of above parameters. The numbers 
are fixed in our simulation, it would not be optimal for all the combinative values oik/ g and 
T/g. But we have tried to used the numbers, which could result in relatively large values of 
P suc and F in the whole range of k/g and T/g. We believe that our numerical investigation 
for the success probability and the fidelity of the two-qubit conditional phase gate subject 
to the spontaneous emission and the cavity decay would be useful for real experiments. 
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Figure captions 

FIG. 1. The level configuration for atom j, where the qubits are encoded in levels \0)j and 
Except that levels and \2)j are coupled resonantly by the cavity mode, all the other 
arrows indicate resonant couplings by different lasers with the corresponding Rabi frequencies 
labeled nearby. 

FIG. 2. Time evolution with the initial state |(|0)i|0) 2 + |0)i|l) 2 + |l)i|0) 2 + |l)i|l)2)|0), where 
k = O.lg and T = O.lg. (a) Pulse shapes and sequences for Qq.i an d ^<r,2- (b) Probability 
amplitudes for states |0) 1 |1) 2 |0), |1)i|<t) 2 |0), |0)i|0) 2 |0), |l)i|0) 2 |0), |l)i|l) 2 |0). (c) Probability 
of all the states containing excited atomic state, i.e., P exc = i^2)i|i) 2 |o) + ^i)i|2) 2 |o) + -^2>i|o) 2 |o> 
, and probability of all the states containing one photon, i.e., P p ht = P\i) 1 \i) 2 \i) + ^i)i|0) 2 |i)- 
(d) Success probability P suc - 

FIG. 3. Dependence of success probability P SU c (upper subplot) and fidelity F (lower subplot) 
on k and T for the initial state 1 0} 1 1 1) 2 1 0} . 

FIG. 4. Dependence of success probability P suc (upper subplot) and fidelity F (lower subplot) 
on k and T for the initial state 1 0> 1 1 0) 2 10) . 

FIG. 5. Dependence of success probability P suc (upper subplot) and fidelity F (lower subplot) 
on k and T for the initial state §(|0>i|0> 2 + |0)i|l) 2 + |l)i|0) 2 + |l)i|l) 2 )|0). 
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